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ABSTRACT 


Compressible  flow  over  a laser  turret  creates  a phase  distortion. 

To  gain  insight  to  this  problem  a model  of  a hemisphei’ical  turret  is 
developed.  The  flow  field  is  obtained  using  a second  order  solution 
based  on  Janzen-Rayleigh  expansion  for  the  compressible  potential  flow 
equation.  Contours  of  constant  phase  distortion  were  calculated;  results 
at  the  critical  Mach  number  are  presented  and  discussed.  It  was  found 
that  the  distortion  at  0°  elevation  angle  was  equivalent  to  a positive 
lens  whereas  at  90°  the  distortion  is  similar  to  a negative  lens.  At 
45°  the  main  effect  is  beam  tilting  into  the  wind.  For  a ratio  of  beam 
radius  to  turret  radius  equal  to  0.5,  the  focal  length  at  beam  elevation 
angle  of  0°  is  16.5km.  At  beam  elevation  of  54°  the  beam  tilt  was  found 
to  be  26.5  microradians . The  results  of  this  paper  are  compared  with 
other  results  reported  in  the  literature. 
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PHASE  DISTORTION  DUE  TO  AIRFLOW  OVER 
A HEMISPHERICAL  LASER  TURRET 

by 

A.  E.  Fuhs*  and  S.  E.  Fuhs** 

I.  INTRODUCTION 

The  index  of  refraction  of  air  is  related  to  mass  density.  Starting 
at  modest  subsonic  Mach  numbers,  e.g.  0.3,  and  extending  to  transonic 
flight,  flew  over  a blunt  laser  turret  causes  significant  changes  in 
density.  The  change  in  density,  which  in  turn  changes  index  of  refraction, 
can  defocus  a laser  beam. 

Variable  index  of  refraction  within  context  of  aeronautical  appli- 
cations is  not  a new  problem;  however,  there  are  two  aspects  which  are 
new.  First,  the  laser  beam  is  large  relative  to  a characteristic  turret 
dimension.  Second,  for  many  airborne  applications  the  laser  telescope 
requires  a turret  which  protrudes  into  the  flow. 

Reconnaissance  aircraft  fly  photographic  missions.  The  compressible 
flow  over  the  aircraft  acts  like  a distorted  lens  and  conceivably  could 
cause  distorted  photographic  images . The  variable  index  of  refraction  is 
not  a serious  problem  for  reconnaissance  aircraft  for  two  reasons.  The 
aperture  of  the  camera  is  small  compared  to  the  characteristic  scale  of 
the  distortions;  in  the  case  of  a laser  telescope  this  in  not  true.  Fur- 
thermore, the  camera  on  a reconnaissance  aircraft  is  not  mounted  in  a 
turret.  To  aim  the  camera  the  pilot  aims  the  aircraft.  Camera  windows 
are  flush  with  the  aircraft  surface. 

Cruise  missiles  with  intercontinental  ranges  may  use  celestial  navi- 
gation which  requires  precise  measurement  of  angular  location  of  stars. 

* Distinguished  Professor  of  Mechanical  Engineering 
**  Student,  California  Institute  of  Technology 
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A distorted  lens  created  by  the  ambient  flow  field  causes  an  error  in 
measurement.  Knowledge  of  the  external  aerodynamics  is  important. 


A model  which  is  useful  for  understanding  the  density  field  in 
the  vicinity  of  a turret  is  the  flew  around  a sphere  or  an  ellipsoid. 

The  plane  containing  an  incoming  fluid  velocity  vector  and  a diameter 
of  the  sphere  does  not  have  flow  components  normal  to  the  plane.  This 
is  a plane  of  symmetry  and  conceptually  can  be  considered  a wall.  The 
flow  over  the  hemisphere  above  the  plane  can  be  calculated  readily. 

For  incompressible  flew,  the  solution  for  the  flow  over  a sphere 

is  well  known.  See  for  example  pages  464-465  of  Milne-Thompson  / 

(2)  (1) 
pages  92-93  of  Lamb  or  pages  339-342  of  Karamcheti.  Techniques 

are  available  to  introduce  the  influence  of  compressibility  or,  stated 

another  way,  the  influence  of  Mach  number.  One  method  is  the  Janzen- 

. . • (4) 

Rayleigh  technique  discussed  briefly  by  Van  Dyke.  Another  brief 

description  is  given  in  Chapter  10  of  Liepmann  and  Puckett/5'*  The 

Janzen-Rayleigf  method  is  discussed  also  on  pages  328-334  of  Oswatitsch/6^ 

(7) 

Lord  Rayleigh  provides  the  second  order  solution  for  a sphere  in  a 
paper  published  in  1916. 
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II.  DISCUSSION  OF  ANALYTICAL  MODELS 


Optical  distortion  due  to  the  external  flow  field  can  be  divided 
into  two  categories.  Viscous  flow  phenomena  fall  into  one  category  and 
include  shear  layers,  laminar  and  turbulent  boundary  layers,  and  shedding 
of  discrete  vortices.  The  other  category  involves  the  inviscid  flow 
field.  This  paper  treats  the  latter  category. 

Laser  turrets  can  be  classified  as  blunt  or  small-perturbation 
turrets.  If  the  turret  has  a surface  with  a normal  vector,  n,  aligned 
parallel,  or  nearly  so, to  the  freestream  velocity  vector,  $ , the  turret 
is  "blunt".  If  the  angle  between  the  normal  vectors  of  the  turret  sur- 
face and  the  freestream  velocity  vector  are  everywhere  large,  e.g.  60° 
to  90°,  then  the  turret  is  "small-perturbation". 

The  turr'et  being  considered  here  is  a hemisphere,  which  is  obviously 
a blunt  turret.  At  the  front  stagnation  point  of  a hemispherical  turret, 

$ • n = - 1.0. 

For  a blunt  turret,  compressibility  effects  set  in  at  a smaller 
freestream  Mach  number,  Mot.  A useful  concept  is  the  critical  Mach  number, 
M^* . At  the  critical  Mach  number,  somewhere  on  the  body  a local  Mach 
number  is  just  sonic.  When  M^  > M^*,  the  flow  becomes  inherently  nonlinear 
and  shock  waves  appear. 

For  a circular  cylinder  with  its  axis  normal  to  the  flow,  the  critical 
Mach  number  is  0.3985.  For  a sphere,  the  critical  Mach  number  is  0.5868. 

A hemisphere-cylinder,  which  is  oriented  with  the  cylinder  axis  parallel 
to  the  freestream  flow  direction,  has  a critical  Mach  number  less  than 

c 

( 8 ) 

0.7.  Hsien  gives  experimentally  determined  pressure  coefficients  for 
a hemisphere-cylinder  over  the  Mach  number  range  0.7  to  1.0.  In  refer- 
ence 9 Hsieh  extends  the  Mach  number  range  from  1.05  to  1.82.  Using 
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from  reference  8,  the  maximum  local  Mach  number  is  1.25;  the  maximum 

local  Mach  number  occurs  at  an  angle  of  approximately  60°  on  the  sphere. 

The  angle  is  measured  from  the  stagnation  point  to  the  location  of 

minimum  C . 

P 

For  small  perturbation  turrets  the  critical  Mach  number  is  larger. 
Calculations  were  made  of  a turret  on  a cylindrical  fuselage/10^  The 
turret  shape  was  determined  by 


R(x,0)  = Rq  + ^-  (1  + cos  (1  + cos  f0) 
R(x,0)  = Rq  for  l < | x | < L and  v/f  < |e|< 


(1) 


where  Rq  is  fuselage  radius;  c,  height  of  turret;  t,  length  of  turret; 
1/f,  the  fraction  of  the  2tt  circumference  occupied  by  the  turret;  and,  L, 
the  separation  between  periodic  turrets.  Table  I summarizes  the  critical 
Mach  number. 


TABLE  I.  Critical  Mach  Numbers  for  Snail-Perturbation  Laser  Turrets 


l/Rn 


L/Rn 


1/f 


e/Rn 


M * 


maximum  slope 
on  turret 


1.005  5.0  .333 


0.05 

.88 

4.5° 

.10 

.82 

8.9° 

.15 

.76 

13.2° 

.20 

.72 

17.4° 

.25 

.68 

21.3° 

.30 

.65 

25.1° 

.35 

.62 

28.7° 

.40 

.59 

32.0° 

.45 

— 

35.1° 

.50 

— 

38.0° 

.55 

— 

40.7° 
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There  are  different  analytical  models  that  can  be  developed  for 
laser  turret  geometries.  Flew  over  blunt  turrets  can  be  determined  using 
expansion  techniques  such  as  the  Janzen-Rayleigh  method  mentioned  ear- 
lier. As  more  terms  in  the  expansion  are  evaluated,  the  more  accurate 
are  the  flew  field  calculations  as  Mach  number  increases.  Upper  bound 
is  the  critical  Mach  number.  For  M > M *,  the  calcinations  will  be 
grossly  in  error. 

There  are  numerical  techniques  which  directly  integrate  appropriate 
differential  equations.  One  method,  which  is  applicable  to  axisymnetric 
flow,  is  the  technique  of  South  and  Jameson. ^ ^ ^ 

For  small-perturbation  laser  turrets  the  equations  can  be  linear- 
ized. Compressibility  effects  in  subsonic  or  supersonic  flow  are  incor- 
porated. Transonic  flew  is  excluded  since  the  transonic  equations  are 
nonlinear  even  for  small  perturbation. 

Figure  1 summarizes  the  regions  in  a body  slope  vs.  Mach  number 
map.  The  flow  within  the  region  ABCD  can  be  described  with  sufficient 
accuracy  using  the  linearized,  small-perturbation  equations.  The  line 
CD  represents  the  critical  Mach  number.  The  line  BC  is  shown  at  a 
slope  of  30°.  The  assumptions  involved  in  the  linearization  of  the  po- 
tential flow  equation  becone  less  and  less  valid  as  the  maximum  slope 
increases.  The  region  CDJ,  which  resembles  an  inverted  triangle,  requires 
solution  of  the  nonlinear  small-perturbation,  transonic,  potential  equa- 
tion. 

The  region  EFGH  defines  the  area  where  Janzen-Rayleigh  technique  is 
most  useful.  The  line  GH, which  defines  the  critical  Mach  number,  is 
shown  dashed;  the  critical  Mach  number  is  a function  of  turret  shape. 

Hence  GH  is  intended  to  qualitatively  suggest  the  upper  bound  for  Mach 
number. 
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III.  JANZEN-RAYLEIGH  EXPANSION  TECHNIQUE 


In  spherical  coordinates,  the  equations  of  motion  are  as  follows: 


2 2 
9u  + v 3u  v 1 9^  _ _ a__  9£ 

9r  r ae  ~ r " ” p 9r  " ~ p 9r 


9v  + v 9v  + uv  _ 1 1 9p  _ a^  9£ 

9r  r 90  r p r 90  pr  90 


(2) 


(3) 


1 


9 , 2 . 

a?  <r  u) 


r sin  O 


9 t 

38  (V 


sin  0)  = 0 


(4) 


Equations  (2)  and  (3)  are  Euler's  equations.  Equation  (4)  is  the  con- 
tinuity equation.  For  flow  over  the  sphere,  there  is  no  dependence  on 
the  azimuthal  angle  in  polar  coordinates.  Consequently,  equations  (2) 

and  (4)  are  two-dimensional  equations  depending  only  on  r and  0. 

• • 2 2 
Multiplying  equation  (2)  by  + u/a  and  equation  (3)  by  - v/a  , 

adding  and  combining  with  equation  (4)  yields 


1 9 


r2 


9r 


(r2u) 


r sin  0 


(v  sin  8)  = K (u2  |H.iil|vtuv|vtv^3vj 
90  az  I 9r  r 90  9r  r 96  j 


The  speed  of  sound,  a,  is  the  local  value  and  changes  from  point  to  point 
in  Hie  flow.  It  is  necessary  to  account  for  variations  in  a.  Using  the 
energy  equation  along  a streamtube,  one  can  show  that 


2 

a 


(1 


u2  + 


) 


(6) 


Inviscid  flew  is  considered.  A potential  can  be  introduced  for  the 
velocity  components  ..  94> 


* ■ ? If  <« 

where  U is  the  frees tream  velocity  at  upstream  infinity.  Introducing 
equations  (7)  and  (8)  into  equations  (5)  and  (6)  yields  equations  in 
terms  of  a potential. 

The  Rayleigh-Janzen  expansion  considers  the  potential  function  to 


be  given  as 


= *0 + M!  *i + m!  t>2 + 


The  solution  for  is  considered  the  first  order  solution.  The  solution 
involving  both  4>q  and  is  called  the  second  order  solution.  In  this 
paper  only  the  second  order  solution  is  considered.  For  ease  of  writing, 
the  following  definition  is  introduced: 

<P1  (10) 

To  digress  briefly  from  the  development  of  the  Janzen-Rayleigh 
expansion,  the  solution  for  the  flow  over  a cylinder  normal  to  the  fluid 

. . . (4)  , 

velocity  will  be  discussed.  Van  Dyke  indicates  the  solution  for  the 

16 

cylinder  has  been  performed  to  Mot  . For  y = 7/5,  the  result  for  maximum 
surface  speed  is 


= 2.00000  + 1.16667M2  + 2.57833M4  + 

(J  OO  00 


+ 25.59041M8  +....+  7450M16 


The  critical  Mach  number  for  a cylinder  is  = 0.39852  + 0.00020. 
In  terms  of  the  potential  function,  the  speed  of  sound  becomes 


S = .i  ♦ F<v  v r’ y> 

d 


(12a) 


where  F is  given  by 


A 


Hr  n-fT^iM  .w-Hi 

ift  ■ 


[■-  (•:•!)] 


(12b) 


With  substitution  of  equations  (9)  and  (10),  the  potential  function 


oecomes 


+ V2*  = M2  } 

oo  7 


9 2 

d)  — o d)  d)  d) 

r Yrr  r2  Wer 


*r  + ?*0  *eej  (13) 


where  all  terms  on  the  right-hand  side  are  determined  from  the  first 

order  potential  function  <j>g.  The  result  of  equating  coefficients  of  like 
2 

powers  in  yields  two  equations 

V2<ju  = 0 (14) 


V2t li  - M2  V24>  - M2  i>2d>  + — p d>  d>.di„ 

oo  T,  oo  Yr  rr  r^  r 0 0r 


?3  *6  V r4  *0  *0©] 


The  right-hand  side  of  equation  (15)  is  known  due  to  the  fact  that 
equation  (14)  can  be  solved  to  yield  . To  emphasize  this  fact,  the 
first  term  on  the  right-hand  side  is  written  explicitly  showing  4>q, 


■©) 


Equation  (15)  is  Poisson's  equation. 
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IV.  FIRST  ORDER  SOLUTION 


The  first  order  solution  is  obtained  from  Milne-Thompson / ^ Lamb/2'* 
(3) 

or  Karamchevi.  It  is 

♦0<r,  0)  = U (r  cos  e + -R32gs--9-  ) (16) 

where  R is  the  sphere  radius.  The  upstream  stagnation  point,  P,  is  at 
P(R,  180°);  and  the  downstream  stagnation  point  is  at  P(R,  0°).  Using 
equation  (16)  one  can  determine  the  derivatives  to  substitute  into  equa- 
tion (15).  These  are 


}>r  = U cos  6 (1  - p-  ) 


<j>  = U cos  6 

rr 


- r U sin  0 


= - r U cos  0 


“ 6 [*  * T £] 

“ t1  ’ £] 


4>q  = - U sin  0 
0r 


The  result  of  combining  equations  (15)  and  (17)  to  (21)  is 

V2^  = |3A(1-X)2  + sin2  0 [4  - | A + 6A2  - ^ A3]  } (22) 

3 3 

where  A = R /r  has  been  introduced. 

Before  proceeding  with  second  order  solution,  a comment  about 
boundary  conditions  is  appropriate.  The  component  of  velocity  normal  to 
the  surface  of  the  sphere  must  be  zero.  At  infinity,  the  velocity  should 
be  U.  Consequently, 


I!i\  = o ■ . o • -1  !M 

Sr  ) r ’ ar  y-  ’ r 11  I. 
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V.  PAYLEIGH  SOLUTION  FOR  SECOND  ORDER  FUNCTION 

(7) 

In  July,  1916,  Lord  Rayleigh  reported  the  solution  for  the 
second  order  function. 

The  Legendre  functions 

P1  = cos  0 (24) 

P3  = |-  cos3  0 - cos  0 (25) 

are  introduced  into  equation  (15)  in  the  following  manner:  First, 

equations  (17)  through  (21)  are  substituted  into  equation  (15);  second, 
terms  are  grouped  so  as  to  form  the  functions  P1  and  P3>  The  result  is 


( 36  R6  9 R9  \ / 

R R3 

24  R6  . 3 

R9  \ 

n 

{'  T¥r+  J Tnrj  pi  + 

6 FT- 

Z~vT  7 

rioj  P3J 

Equation  (26)  is  in  the  form 


72<j>1  = z z r'  PC 


where  C^  are  zero  except  for  m = 7 and  10;  all  C2m  are  zero;  and  C3m 
are  zero  except  for  m = 4,  7,  and  10.  Non- zero  C are  the  appropriate 
constants  appearing  in  equation  (26).  The  Legendre  differential  equation 
is  dp 

sht  m [sin  8 ar]  * " (ntl>  pn  = 0 (28> 


Define  the  function  <p  which  satisfies  the  equation 

nm  n 


a2<f>  „ 3<J>  . 

r_rnP  C = V2*  = — ™ + 1 _™L  + 1_ 

n nm  nm  3r2  r 3r  r2sin 


1 3 f.  „ 

sin  0 30  S;Ln  6 


p = r2_m  P A 
nm  n nm 


(30) 


Consequently, 


(31) 


9<f) 


IE  = (2-m)  r1-m  P A 
3r  n nm 


and 


924> 


-4^-  = (2-m)  (1-m)  r'm  P A 
9rz  n nm 


Combining  equations  (29)  and  (32)  gives 


C = (2-m)  (1-m)  r_m  P A 


n nm 


n nm 


— m+2 

+ 2 (2-m)  r-mP  A + — . — - 2- 
n nm  r2  sin  0 96 


[3P  i 

sin  6 SfT_ 


nm 


Substitution  of  equation  (28)  into  equation  (33)  yields 

’n  £(2-m)  (1-m)  + 2 (2-m)  - n (n+1)^ 


r m P C = A r m P 
n nm  nm 


As  a result 


nm 


nm  (2-m)  (3-m)  - n (n+1) 


Using  equations  (26),  (30),  and  (35),  that  part  of  the  second 
function  identified  as  the  particular  solution  is  obtained 


R6  . R9  1 p 

3R3 

3R6  3R9  1 ) 

“ 5FT  + TO'  + p; 

3 " lDr?  " 

I0r5  r75r®"  j 

The  homogeneous  solution  is 


!>1  = Ar  pi  + BP1/r2  + Cr  P3  + D P3/r4 


(32) 


(33) 

(34) 

(35) 
order 

(36 

(37) 


At  infinity  the  potential  must  be  finite  and  the  velocity  must  be  U; 
the  function  ^ of  equation  (16)  yields  U.  Consequently  both  A and  C 
equal  zero.  At  the  surface  of  the  sphere 


r 


I 


(38) 


From  the  requirements  of  equation  (38)  combined  with  equations  (36) 


and  (37) 
and 


(39) 


D = 


_ 27R5U 


"55" 


M2 


(40) 


The  final  expression  for  the  potential  function  is 


♦ = VM®* 


Hu  [r ' 


3R3  27R5  3R6  , 3R9 

10r2  55r4  ~ 10r5  176r8 


♦ .51 


®^3r2  ~ 5r5  24r8|  *1 


'] 


]’•( 


(41) 


The  first  term  in  brackets  is  4>Q ; and  the  second,  long  tern  is  the  second 


order  contribution. 


VI.  DEVELOPMENT  OF  PHASE  DISTORTION 


Use  of  Potential  Function  to  Obtain  Index  of  Refraction 
The  index  of  refraction  and  density  are  related  by 


n = 1 + k*  = = 1 + < (4 

poo  PSL  p® 

• —4 

where  tc'  is  a constant  with  value  of  about  2 x 10  . Note 

k = < ' p^/ pc^.  Tne  flow  is  assumed  to  be  isentropic.  Consequently 


1 + ^ 

p _ 2 °° 

p®  1 + M2 


The  local  Mach  number  is  given  by 


M2  = p[(i)2^  n)2] 


In  equation  (42)  p is  the  local  density,  p^  is  the  freestream  density, 
and  pC7  is  the  sea  level  density  on  a standard  day.  By  combining  equa- 
tions  (6),  (12a),  (41),  (42),  (43),  and  (44),  an  expression  is  obtained 
for  the  index  of  refraction. 

The  turret  has  a characteristic  length;  for  the  sphere  this  is 
the  radius  R.  The  laser  beam  has  a characteristic  size  which  is  the 
diameter  D or  radius  R' . The  wavelength  of  the  laser  radiation,  X, 
is  also  a characteristic  length  as  will  be  demonstrated  in  subsequent 
discussion.  Two  nondimensional  lengths  can  be  formed.  Since  R occurs 
prominently  in  the  equations  for  the  potential,  R is  chosen  as  the 
reference  length.  With  that  choice,  X/R  and  D/R  are  two  of  the  para- 
meters specifying  the  problem. 
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Phase  distortion  is  a function  of  the  direction  that  the  beam  is 
pointed.  Define  e and  a as  elevation  anc  azimuth  angles.  For  the 
hemispherical  turret,  symmetry  exists  about  the  diameter  of  the  sphere 
which  is  parallel  to  the  freestream  velocity  vector.  In  general,  two 
angles,  e and  a,  are  needed  to  define  the  beam  direction;  however, 
for  the  hemispherical  turret  only  one  angle  is  needed.  Looking  ahead, 
the  angle  0'  occurring  in  equations  (45)  and  (46)  will  be  used  as  the 
beam  angle;  see  also  Figure  2. 

Geometrical  Relations  for  Calculation  of  Phase  Distortion 


Figure  2 shows  the  geometry  of  the  problem.  The  radius  of  the 
sphere  is  R.  The  beam  is  elevated  at  an  angle  0 ' . The  z axis  is  the 
same  as  the  polar  axis  for  polar  coordinates. 

To  locate  points  within  the  beam,  three  variables  are  used;  s is 
distance  along  the  beam,  as  shown  in  Figure  3;  a is  an  angle  measured 
clockwise  from  the  windward  edge  of  beam;  and  R'  is  the  radial  distance 
from  the  beam  centerline  or  axis. 

The  flow  properties,  i.e.,  density  and  local  Mach  number,  are 
filiations  of  r and  0.  These  are  polar  coordinates.  Consequently,  it  is 
necessary  to  express  r and  0 in  terms  of  0',  s,  a,  and  R.  The  equations 

^ 2 T . I2  2 2 

r = ! (R  + s)  sin  0*  - R'  cos  a cos  0’J  + R'  sin  a 

r -i  2 (45) 

+ I (R  + s)  cos  0'  + R’  cos  a sin  6’ 


0 = tan 


^ ^j(R  + s)  sin  0'  - R*  cos  a cos  0'j  2 + R'2  sin  2a 
|^(R  + s)  cos  0*  + R*  cos  a sin  0'j 


Equations  (45)  and  (46)  can  be  used  to  evaluate  equations  (41)  to  (44). 
The  index  of  refraction  is  given  by  equation  (42).  Define  N as  the 


16 


number  of  waves  along  a distance  s.  Hence 


dN  = (47) 

where  A is  the  local  wavelength.  The  local  wavelength  is  a function  of 
0' , a,  R' , and  s.  It  is 


v c 

I~  HF 


As  a result,  the  number  of  waves  is 


A (1  + k) 
00 

(l  + Kp/p„ 


A u + <5 


(1  + k — )ds 
Pm 


A reference  is  needed.  The  reference  which  is  convenient  is  the  number 
of  waves  along  the  beam  axis;  designate  that  quantity  as  Nq.  Using  this 
reference 


N-No  = r 


Cl+K  / (1+k  2-)ds  - J (l+<^)ds  - j (1+Adsj  (50) 


fJote  the  density  on  the  beam  centerline  for  -s'<_  s _<  0 is  denoted  by  pr, 
a reference  density.  For  0 <_  s < <*>,  the  density  on  the  centerline  is 
Pq  which  is  a function  of  s.  In  equation  (50)  the  value  of  unity  in 
each  integrand  can  be  cancelled  giving 


The  second  integral  in  equation  (51)  has  been  termed  the  GAP  INTEGRAL 
since  it  enters  the  equation  as  a result  of  the  gap  between  s = -s'  and 
s = 0.  See  Figure  3.  In  the  calculatlors  which  were  conducted,  the  gap 
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integral  was  evaluated  by  assuming  = pm;  this  is  an  arbitrary  assump- 
tion. Note  that  in  equation  (51)  wavelength  is  nondimensionalized  by  R 
to  form  WR.  The  distance  along  a ray  is  also  nondimens ional,  i.e., 
ds/R. 

A computer  program  has  been  developed  for  the  HP9830  which  calculates 
N - Nq  as  a function  of  sphere  radius,  R;  elevation  angle,  0';  radius 
within  the  beam,  R' ; angle  within  the  beam,  a;  freestream  Mach  number, 

Mi  wavelength,  A^;  constant  in  density  equation  (42),  ratio  of  heat 
capacities  for  air,  y;  and  freestream  speed  of  sound,  am.  Note  that 
it  is  not  necessary  to  q^ecify  p^  since  all  calculations  depend  on  den- 
sity ratio;  see  equations  (48)  and  (51).  Values  used  in  the  calculations 
are  as  follows: 

R = 0.9144  meter 

_6 

A =3.8  micron  = 3.8  x 10-  meter 
00 

k'  = 2 x 10~4 
a = 342  meter/ sec 
y = 1.4 

M , R' , a,  and  91  have  been  varied.  Calculations  have  been  performed 
to  the  extent  necessary  to  plot  phase  distortion  maps. 
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VII.  RESULTS  OF  CALCULATIONS 
Graphical  Presentation  of  Results 


TABLE  II.  MAXIMUM  VALUE  OF  GAP  INTEGRAL 

beam  elevation  angle  maximum  value  of  gap  integral 

in  deciwave  lengths  and  location 


0' 

within  beam 
(gap  integral) 

a 

0° 

+ 36.8 

ft 

o 

00 

rH 

+ 65.4 

0° 

36° 

+ 67.8 

0° 

O 

-If 

lO 

- 83.6 

180° 

72° 

- 86.3 

120° 

O 

o 

CD 

- 86.6 

90° 

; 


Evaluation  of  Phase  Distortion  as  a Function  of  Distance  Along  the  Beam 
An  important  question  is  the  rate  of  decay  of  the  integrand  in 
equation  (51).  Calculations  were  made  of 

3(N  - NQ) 

?(s/R) 

which  is  the  phase  distortion  per  unit  distance  along  the  beam.  The 
phase  distortion  per  unit  distance  is  a function  of  distance  along  the 
beam.  Results  of  the  calculation  are  shown  in  Figure  10.  The  phase 
distortion  is  the  area  enclosed  by  one  of  the  curves.  Most  of  the  phase 
distortion  occurs  within  a length  equal  to  one  turret  radius. 
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Another  parameter  of  considerable  interest  is  the  dependence  of 
phase  distortion  on  Mach  number.  Qualitatively,  isodensity  contours  at 


M^  < M*  have  an  appearance  the  same  as  for  Moo=  M*.  Quantitatively,  there 
are  changes  in  both  the  shape  of  an  isocontour  of  constant  phase  shift 
and  the  magnitude  of  the  phase  shift.  Table  III  compares  magnitudes 
of  distinctive  features  of  the  naps  for  two  Mach  numbers.  Table  III 
should  be  studied  simultaneously  with  the  appropriate  figure  for  6 1 . 

Phase  distortion  curves  based  solely  on  a first  order  calculation, 
i.e.,  using  4>q  in  equation  (9)  only,  should  have  the  same  shape  for  all 
M^.  A calculation  of  phase  distortion  based  on  a second  order  solution, 
i.e.,  using  both  <J>q  and  <J>^  in  equation  (9),  should  be  a function  of  M^. 
Figure  11  shows  the  change  in  shape  of  a phase  distortion  curve  with 
Mach  number.  The  curves  are  drawn  for  6'  = 18°:  M , either  0.300  or 
0.587;  R'/R  from  0 to  0.5  and  for  a = 0°  and  a = 180°.  The  phase  dis- 
tortion has  been  normalized  with  the  value  at  R'/R  = 0.5  and  a = 180°. 
Differences  in  the  curves  are  due  to  Mach  number  effects. 
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0’ 

feature  of  isocontour  map 

value  at  M = 0.300 
00 

M = 0.587 
00 

0° 

maximum  <J>  at  beam  edge 

-41.5 

-160 

00 

o 

maximum  <p  at  positive  peak 

+15 

+ 60 

minimum  4> 

-69.6 

-290 

36° 

maximum  $ 

+44.2 

+240 

minimum  <(> 

-80.3 

-345 

54° 

maximum  <{> 

+85.2 

+370 

minimum  <p 

-54.4 

-268 

72° 

neximum  <p 

+83.7 

+386 

minimum  4)  at  negative  valley 

-14.1 

- 74.1 

90° 

maximum  <t>  at  a = 0° 

+49.1 

+ 223 

M.«  0.300 


NONDIMENSIONAL  RADIUS 

Figure  11.  Fhase  Distortion  for  Two  Mach  Numbers. 
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Interpretation  of  Results 

Phase  shift  for  6'  = 0 as  shown  in  Figure  4 consists  of  a series 
of  concentric  circles.  Since  the  phase  shift  is  negative,  the  beam 
is  being  focused.  Concentric  curves  with  negative  phase  distortion 
gives  a positive  focus. 

The  contours  of  constant  phase  shift  for  0'  = 54°  are  almost 
straight  lines.  The  algebraic  signs  of  phase  shift  contours  indicate 
the  beam  is  being  tilted  in  a direction  opposite  to  the  relative  wind. 
The  angular  tilt  is  given  roughly  by 

(tilt  angle)  = ^N-NQ)0O  - (N-N0>180oj  ^ (52) 

Inserting  the  appropriate  values  (values  for  N-NQ  were  obtained  from 
Table  III)  gives 


(tilt  angle)  = 


(370  + 268)  (3.8  x 10-6m) 
(2)  (.4572m)  (100) 


26 .5  microradians 


Tilt  of  a beam  does  not  cause  a serious  problem. 

These  results  for  tilt  angle  can  be  compared  with  similar  work  by 


Cook. 


(12) 


(12) 


Captain  Cook  used  the  following  values  in  a sample  calculation 
(notation  has  been  adjusted  to  match  this  report) : 


R'  = 0.25m  A =10.6  micron 

oo 

R = 0.71m  k'/P"  = 0.22  cm3/gm 

M = 0.5  k'  = 2.7  x 10~4 

00 

For  0'  = 54°,  a beam  tilt  of  63  microradians  was  calculated.  When  adjust 
ments  are  made  for  different  parameters  used  to  calculate  the  two  results 
i.e.,  63  and  26.5  microradians,  the  answers  agree  quite  favorably. 

It  is  possible  to  determine  an  approximate  focal  length,  F . The 
subscript,  ct,  indicates  the  focal  length  in  a plane  defined  by  a of 
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Figure  2.  The  appropriate  equation  is 


r-  * TfcTi  <n-no>r,  - 2(m0>r, 


Equation  (53)  was  evaluated  for  a = 0 and  for  the  following  values  for 
other  parameters 

p = 1.23  x Id-3  gm/cm3 


X = 3.8  x 10'u  m 


M = 0.587 


AR'  = 0.1524  m 


R = 0.9144  m 


(N-Nq)  at  =0.799  wavelengths 
(N-Nq)  at  R|  =0.214  wavelengths 
0'  = 0° 

A focal  length  of  16.5  km  was  found. 

Before  proceeding,  some  additional  information  concerning  equation 
(53)  is  appropriate.  Equation  (53)  is  a specialization  of  the  equation 

1 _ 32I 


where 


Fx-  557 


to 

"I  " 


For  the  problem  at  hand  equations  (54)  and  (55)  yield 


1 K 

Fx  ‘ TtoT2' 


p2~p0 


po-po 


where  p2  = p(x2) , p^  = p(x1),  etc.  It  is  obvious  that  the  third  integral 

(12) 

is  zero.  Equation  (56)  leads  to  equation  (53).  From  Cook's  report  for 
the  previously  cited  values  for  R,  R' , M^,  etc.,  a focal  length  of  about 
6 km  is  calculated.  Adjusting  the  value  for  values  of  R,  R',  M , etc., 
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used  in  the  calculations  of  this  report,  good  qualitative  agreement 
is  obtained. 

(13) 

Wolters  and  Laffay  have  calculated  focal  length  for  an 
aperture  of  a turret  on  a F-15.  The  turret  geometry  was  a cylinder 
with  axis  parallel  to  the  freestream;  radius  of  the  cylinder  was  0.5  m. 
When  viewed  from  a position  normal  to  the  axis,  a cylinder  appears  as  a 
rectangle.  The  square  corners  were  rounded  with  a radius  of  0.3  m. 

Without  rounding  of  the  comers,  the  cylinder  viewed  from  its  end  would 
appear  as  a disc  with  radius  0.5m.  With  the  rounding  of  the  comers , 
the  disc  is  reduced  to  0.2  m;  the  disc  is  normal  to  the  freestream 
velocity  vector. 

(13) 

Wolters  and  Laffay  determined  the  focal  length  of  the  lens 
created  by  density  variations  within  the  external  flow.  The  following 
parameters  were  used  in  their  calculations. 

R'  = 0.17  m X = 10.6  x 10"6  m 

oo 

R„  = 0-5  m K'/p  = 0.223  x 10~3  m3/kg 

M = 0.587  0'  = 0 

OO  w 

altitude:  20,000  ft. 

For  these  conditions  a focal  length  of  18.2  km  was  determined.  Since 
the  distortion  varies  linearly  with  the  ambient  density,  a fact  which 
follows  from  equation  (42),  the  focal  length  can  be  adjusted  to  sea  level 
by  an  adjustment  to  density.  For  the  same  conditions  as  above  except  at 
sea  level,  the  focal  length  would  be  8.2  km.  A value  8.2  km  is  consistent 
with  results  obtained  in  this  report,  16.5  km. 

Having  compared  in  a rough  way  calculated  results  from  references 
(12)  and  (13)  with  calculated  results  from  this  paper,  the  relative 
advantages  and  disadvantages  of  each  model  will  now  be  compared.  Reference 
(12)  used  3l/3x  and  3l/3y  to  predict  tilt;  32I/3x2  [see  equation  (54)]  and 
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32I/3y2  were  used  to  predict  focal  length.  These  were  evaluated  at  the 
beam  center.  An  analytical  solution  was  obtained.  Reference  (12)  also 
used  a first  order  solution  for  the  flow  field,  which  is  equation  (16). 
This  paper  uses  a second  order  solution  as  given  by  equation  (41).  Hence 
a more  accurate  flow  field  is  used  in  this  paper.  Reference  (12)  yields 
a solution  in  simple  form;  this  paper  requires  a desk  top  computer,  eg. 
HP9830,  or  else  a high  degree  of  patience  if  a HP55  were  to  be  used. 

In  contrast  to  reference  (13)  which  used  a numerical  solution  to 
the  flew  field,  the  second  order  analytical  solution  gives  the  flow  field 
to  any  degree  of  spatial  resolution  required.  Numerical  solutions  give 
results  in  discrete  steps.  Based  on  Figure  5 of  reference  (13)  the  ratio 
of  step  size  to  cylinder  radius  is  0.22.  Figure  10  of  this  paper  shows 
that  approximately  one-half  of  the  phase  distortion  is  generated  in  a 
step  size  0.22.  Of  course  a smaller  step  size  can  be  used.  Numerical 
techniques  have  the  advantages  of  providing  solutions  for  complex  turret 
geometry.  This  paper  is  applicable  only  to  a hemispherical  turret. 

Look  at  Figure  9 for  which  6'  = 90°.  The  contours  have  a symmetry 
about  the  a = 0°--180°  line  and  about  the  a = 90°--270°  line.  The  phase 
distortion  has  positive  values  which  give  a negative  lens.  The  focal 
length  for  a - 0°  is  considerably  less  than  for  a = 90°.  At  9'  = 90°, 
the  beam  is  defocused. 


It  is  apparent  from  Figures  4 to  9 that  the  phase  distortion  is  a 
strong  function  of  elevation  angle  6 ' . Compensation  for  the  phase  dis- 
tortion is  possible,  at  least  partially.  However,  the  shape  of  the  pri- 
mary telescope  mirror  required  for  compensation  is  a function  of  both 
elevation  angle  and  Mach  number.  Such  compensation  should  fall  within 
the  capability  of  adaptive  optics  techniques. 
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Both 


To  discuss  the  altitude  dependence,  recall  k = k'p^/o^. 
k'  and  psl  are  constants  for  a given  wavelength  and  specified  atmospheric 
conditions.  Consequently,  k varies  directly  with  pm.  Now  examine  equa- 
tion (51).  Since  k <<  1,  the  term  1 + k is  essentially  unity;  it  follows 
that  N-Nq  varies  directly  as  k.  Phase  distortion  changes  as  pw.  Obviously 
as  altitude  increases,  phase  distortion  decreases. 
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VIII.  CONCLUDING  REMARKS 


Compressible  flow  over  a sphere  has  been  obtained  forming  the  basis 
for  a model  of  a blunt  hemispherical  laser  turret.  Phase  distortion  has 
been  calculated  for  several  frees tream  Mach  numbers.  Phase  distortion 
data  are  presented  for  M^  equal  to  the  critical  Mach  number.  Calculations 
are  consistent  with  those  of  references  (12)  and  (13). 

The  phase  distortion  is  a function  of  both  elevation  angle  and  Mach 
number.  Hence  optical  compensation  requires  a variable  shape  primary 
mirror. 

The  model  of  a hemispherical  laser  turret  provides  insight  to  the 
behavior  of  phase  distortion  as  a function  of  Mach  number.  See  Figure 
11.  Also  refer  to  Table  III. 

Another  valuable  aspect  of  the  calculations  is  the  decay  of  den- 
sity perturbation  with  distance  along  the  beam.  See  Figure  10.  It  is 
apparent  that  the  distortion  is  fairly  localized.  Within  a distance  of 
s/R  = 1.0,  over  99%  of  the  distortion  occurs. 

The  formulation  clearly  shows  the  altitude  dependence  of  phase 
distortion. 
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